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Abstract. The structure, lattice dynamics and mechanical properties of the 
magnesium hydroxide have been investigated with static density functional theory 
calculations as well as ab initio molecular dynamics. The hypothesis of a 
superstructure existing in the lattice formed by the hydrogen atoms has been 
tested. The elastic constants of the material have been calculated with static 
deformations approach and are in fair agreement with the experimental data. 
The hydrogen subsystem structure exhibits signs of disordered behaviour while 
maintaining correlations between angular positions of neighbouring atoms. We 
establish that the essential angular correlations between hydrogen positions are 
maintained to the temperature of at least 150 K and show that they are well 
described by a physically motivated probabilistic model. The rotational degree of 
'"^ ' freedom appears to be decoupled from the lattice directions above 30K. 
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1. Introduction 

The magnesiuiTL hydroxide — Mg(0H)2 — is a very simple mineral, called brucite, 
often regarded as a good example of a material with structurally bound OH group. 
The large number of known silicate phases with such structurally bound OH units 
[TJ 12] motivates high interest in the role played by the OH groups and the physical 
properties induced by them of this type of crystals in the Earth interior. Brucite was 
a subject of numerous experimental H [3 HI 13 H [8] and theoretical [9l [TOl [TTl [12] 
H \ studies in the past. Despite its simple structure, its physical properties remain still 

unexplored — in fact even its crystallographic symmetry has not been fully established 
[13j , although the question about its structure has been addressed several times in the 

past [iiiniiniiii]. 

The crystal lattice of the magnesium hydroxide is a simple, layered, hexagonal 
structure with layers of magnesium interlaced with layers of hydroxyl groups. This 
basic structure was firmly established long time ago [14, as having PSml symmetry 
with 0-H bond aligned on a three-fold axis (see figure [T]). Further research [15] showed 
unexpectedly high thermal displacements of hydrogen atoms. The models proposed 
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Figure 1. Schematic view of the structure of a brucite crystal. The figure shows 
2x2x2 supercell used in the numerical calculations, while the primitive unit cell 
is indicated by the (smaller) grey box. Dark grey (red) and light grey (green) balls 
show the positions of O and Mg atoms, respectively. The labels of the hydrogen 
atoms shown by light (white) balls are used to analyse their structural correlations 
in the text. 



until now [51 [71 [T31 [THl HZl [H] deal with temporary and spatially averaged positions 
of the hydrogen subsystem and suggest that the hydrogen atoms are in fact displaced 
from the high symmetry points on the three-fold axis to form a new lattice with 
P3 symmetry and averaged 1/3 occupancy of three equivalent 6z Wyckoff positions 
(x,22;,z), with either x > 1/3 (so called XGT arrangement) [9l H [IS [13 [18], or 
X < 1/3 (so called XLT arrangement) as proposed by Megaw [TB]. The schematic 
placement of atoms in this larger unit cell is presented in figure [2l The case for the 
XGT arrangement is further strengthened by the neutron diffraction experiment of 
Desgranges et al. [12] indicating that protons are displaced into XGT positions even 
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Figure 2. Schematic atomic positions in the P3 structure of a brucite crystal. 
The primitive unit cell is the size of v^ X v^ X 1 of the unit cell of the P3m 
structure shown in figure [l] The corners are occupied by Mg atoms (grey, red 
balls), and each H atom (dark, blue ball) is accompanied by an O atom (light, 
green ball). The labels of the hydrogen atoms are used later in the text to 
investigate their spatial correlations. The angles shown here for h5 atom represent 
its relative angular coordinate tp with respect to the reference atom HI, and an 
absolute angular coordinate 5. 



at ambient pressure. 

All the above models could reflect true structure of the brucite crystal but 
due to the performed spatial and temporal averaging of the structure they may 
miss the dynamical aspect. In the present paper we deal with this aspect of the 
magnesium hydroxide structure. Firstly, we are going to show results of the static, 
elastic constants calculations for the crystal followed by the results of the density 
functional theory (DFT) molecular dynamics (MD) calculations of the brucite crystal's 
structural properties. Secondly, we analyse the dynamical aspects of the structure and 
demonstrate that the hydrogen atoms are not frozen in their equilibrium positions but 
fluctuate strongly, displaying some unexpected dynamical structure characterized by 
long range spatial correlations. 

The paper is organized as follows. In section[2]the calculation method is presented. 
It serves to derive elastic constants of the magnesium hydroxide which are reported 
and compared with available experimental values in section |31 Next, in section HTTl we 
investigate the dynamical model of the hydrogen structure, showing that the short- 
range correlations persist and relative hydrogen positions may be described by a simple 
probabilistic model presented in section |42|) . The paper is concluded in section [5l 
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Table 1. Elastic constants for increasing pressure P (all quantities in GPa) 
obtained for the high symmetry (P3ml) unit cell of the magnesium hydroxide. 



p 


Cii 


C33 


C12 


Cl3 


C44 


1.4 
2.1 

2.7 


126.6 
131.2 
135.5 


534.0 
541.8 
549.5 


112.7 
116.1 
118.9 


39.6 
40.8 
41.5 


61.8 
64.3 
66.9 



2. Calculation setup 

We have performed all calculations of the structure and elastic properties of the 
magnesium hydroxide using VASP,[5n] a DFT based code employing the generalized 
gradient approximation (GGA) ^T with Perdew, Burke, and Ernzerhof (PBE) [22] 
projector augmented- wave (PAW) pseudopotentials [23l [21]. The reciprocal space 
integration has been carried out over the set of special points generated according to 
the Monkhorst-Pack scheme [24] . 

The static calculations have been carried out with a single hexagonal unit cell 
shown by a grey box in figure [TJ a v3 x v3 x 1 supercell (figure [5]) and reciprocal 
space sampling based on 2 x 2 x 2 lattice. The calculation of elastic constants has used 
the standard finite deformation method described in detail in the previous publications 
[2S1[2S]. We have applied displacements between {0.1%; 0.1°} and {1%; 1°} for the 
lattice vector lengths and angles, respectively. 

The schematic structure of the crystal in the 2x2x2 supercell is presented in 
figure [1] together with the adopted numbering scheme of hydrogen atoms used later in 
section U in the analysis of inter-atomic correlations. The MD calculations have been 
conducted with both -\/3 x \/3 x 1 and 2x2x2 unit cells with the standard energy cut- 
off, Ecut=400eV, determined by the used pseudopotentials. We have performed several 
3 ps runs with 1-1.5 fs time step for temperatures ranging from 10 K to 150 K in micro- 
canonical ensemble (constant total energy) with at least 10 ps thermalization runs 
carried out with Nose thermostat technique [27]. The trajectories collected from the 
MD runs were further analysed to extract positional correlations using the computer 
code developed for this purpose. 

To check the influence of the supercell size on obtained results we have calculated 
— using standard direct method [551 HH] — the force constants matrix and full lattice 
dynamics for the larger (3x3x3) system as well as the smaller one used in the 
MD calculations (2x2x2). The obtained exponential dependence of force constants 
on distance shows that the inclusion of the next layer of unit cells in the system 
contributes less than 1/1000 of the overall interaction between atoms in the system. 
Furthermore, the dispersion curves and phonon density of states spectra, calculated 
for both (2x2x2 and 3x3x3) systems, display no significant difference, which 
means that the inclusion of the additional layer of unit cells does not modify the lattice 
dynamics of the system in any significant way and the conclusions concerning lattice 
vibrations drawn from the present calculations are robust. 

3. Elastic constants 

We have completed a series of static, elastic constants calculations using various 
sizes and symmetries of the unit cells and atomic positions: high-symmetry (P3ml) 
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Table 2. Elastic constants of the magnesium hydroxide for increasing pressure 
P (all quantities in GPa), measured by Jiang et al. 1301 with Brillouin scattering. 



p 


Cii 


C33 


C12 


Cl3 


C44 


10-* 


159.0(16) 


49.5(7) 


43.3(17) 


11.1(25) 


22.8(4) 


1.1(1) 


164.1(10) 


59.2(7) 


45.5(8) 


12.3(8) 


25.2(3) 


2.1(1) 


169.3(9) 


73.2(7) 


46.2(8) 


16.6(7) 


28.9(3) 


3.5(1) 


179.6(10) 


94.0(8) 


52.0(9) 


22.1(8) 


33.4(3) 


4.8(1) 


184.3(8) 


110.0(6) 


53.1(6) 


28.2(5) 


36.8(3) 



Table 3. Elastic constants for increasing pressure P (all quantities in GPa) 
obtained for the low symmetry (P3) unit cell of the magnesium hydroxide. 



p 


Cii 


C33 


C12 


Cl3 


C44 


0.01 


130.6 


48.5 


70.3 


10.0 


20.4 


0.96 


133.5 


51.9 


67.8 


9.8 


24.2 


2.0 


143.4 


82.6 


73.4 


17.7 


29.3 


3.0 


150.0 


93.5 


76.3 


15.7 


32.9 


4.0 


153.9 


104.2 


76.1 


22.1 


36.6 



experimental unit cell [14], low-symmetry (P3) y/3 x y/3 x 1 supercell, both with 
standard (1%) and small deformations (0.1%). In general, some results presented in 
table [1] for the high-symmetry unit cell show fairly good agreement with experimental 
data |30] , but others exhibit certain peculiarities which motivated our further 
investigation. While the obtained values of Cn are rather close to experimental ones, 
the calculated elastic constants C12, C13 and C44 are systematically about twice higher 
than the experimental values (table [2]). The largest difference close to one order of 
magnitude was obtained for C33 at P ~ 1 GPa. 

During calculations we have observed that the system shows large sensitivity 
of elastic constants to external stress applied to the unit cell. In some cases non- 
hydrostatic stress of approximately 1-2 GPa induced twofold change in C33 and C13 
constants. 

The discrepancy between theoretical and experimental values of elastic constants 
can be partially resolved by careful relaxation of the internal degrees of freedom of 
the system in the deformed state. The standard procedure calls for cell deformation 
while keeping the same fractional positions of atoms followed by the internal degrees of 
freedom relaxation. Since the hydrogen atoms start at high symmetry positions on the 
three-fold axis the straight-forward minimization procedure keeps this symmetry. Only 
after lifting this symmetry constraint, the true ground state obtained after deformation 
(with lower total energy) could be reached. The elastic constants obtained following 
this procedure become fairly close to the results of Brillouin scattering experiment 
summarized in table [51 As a result of this relaxation the hydrogen atoms are displaced 
from their starting high-symmetry positions {^,^, z) and (|, ;j, 1 — z). 

Prompted by the above anomalies and concepts presented in the published 
research on the structure of Mg(0II)2, we have calculated the same properties for the 
larger \/3 x a/3 x 1 unit cell with P3 symmetry and hydrogen atoms displaced into XGT 
positions as suggested recently by Mookhejee and Stixrude [T^ . This suggestion agrees 
also with the earlier results in the literature [H [3 [ITj [18] . Indeed, the values obtained 
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Figure 3. Absolute angular hydrogen distribution in P3 unit eell. Top panels 
show results obtained from MD run at 244 K, and bottom ones at 46 K. The 
maxima at angle ip ~ 120° obtained for H2 (H5) atoms are 2.05 and 1.2 ( 1.2 and 
1.57 ) at high and low temperature, respectively. Positions of hydrogen atoms 
with labels H1,H2,- ■ -,116 are shown in figure 1. 



for this phase and presented in table [3] show much better agreement with experiment 
(table [2|) . We have also observed considerably weaker sensitivity of calculated values 
to the external stress applied to the system. 

Moreover, during optimization of the magnesium hydroxide structure, we have 
observed that positions of hydrogen atoms around the three-fold axis are not very 
well constrained in the system, i.e., they probably move in a very shallow potential 
well. This observation is in concordance with previous experimental evidence |15j . 
It suggests that the structure of the magnesium hydroxide crystal may be more 
complicated than it is currently accepted. 



4. Dynamical structure 

4-1. Correlations between hydrogen positions 

In order to explore further the issue of dynamical hydrogen structure in the magnesium 
hydroxide, we have used the DFT-based MD technique. The first calculation setup 
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Figure 4. Relative angular hydrogen distributionin P3 unit cell. The Hi box 
shows absolute angular distribution of the Hi atom. Top panels show results 
obtained from the MD run at 244 K, and bottom ones at 46 K. Positions of 
hydrogen atoms with labels H1,H2,- ■ -jHG are shown in figure 1. 



used the -s/S x %/3 x 1, P3 unit cell proposed by Mookherjee and Stixrude ^13,. The 
data have been collected from measurement runs of 3 ps length each preceded by at 
least 10 ps thermalization run, for several temperatures between 15 K and 250 K. 

The examination of the MD trajectories shows large amplitudes of hydrogen 
movements even for runs performed at low temperature. This result corresponds to 
large thermal displacements observed in the experiment J15j . On the one hand, some 
experimental results |19| . showing clear symmetry of the hydrogen locations, may be 
seen as contradicting such observation, while on the other hand, correlations between 
the hydrogen movements could explain this discrepancy. 

To resolve this dilemma we have developed a method to quantify angular 
correlations of hydrogen atoms. We extracted the correlation between positions of 
hydrogen atoms with respect to each other as well as with respect to the fixed reference 
crystal lattice. The results are presented in figure [3] for absolute hydrogen angular 
distribution and in figure |4] for relative hydrogen angular correlations. All these runs 
have been performed in the v^ x v^ x 1 unit cell. The angle (p is an angular position 
of the hydrogen atom around the three-fold axis, measured starting either from the 
direction of the crystal vector A for absolute angular distributions (red histograms). 
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Figure 5. Absolute angular distribution of hydrogen atoms in the 2x2x2 
supercell. The results are obtained from a 6 ps MD run at 16 K. Positions of 
hydrogen atoms with labels H1,H2,- ■ -,1116 are shown in figure 1. 



or from the direction selected by an arbitrary reference hydrogen atom (Hi [i i zi] in 
our case) for relative angular distributions (blue histograms). 

The high peaks seen in figure [3] suggest that there may be a lock- in of H2 and H5 
hydrogen atom positions with respect to the crystal lattice. The lattice lock-in appears 
to be still present at 244 K (upper panel in figure [3|) but we have to mention fairly 
strong and uniform background of other positions in the histogram. In fact, the most 
of the area in these histograms is outside the peaks. Thus, hydrogen atoms do not 
spend most of their time in the vicinity of their high symmetry positions. This means 
that, despite distinct peaks visible in the histograms, the angular positions of hydrogen 
atoms are not locked with respect to the crystal lattice orientation. Some of our results 
obtained for even lower temperatures (T « 10 K) indicate that the true lock-in with 
the crystal lattice may appear at or below T = 15 K — at this temperature we have 
observed Gaussian-shaped peaks and diminished background in the histograms (not 
shown). Unfortunately, due to the difficulties of the MD calculations performed at 
low temperature we were unable to collect enough evidence to confirm this indication. 

The relative angular distribution results displayed in figure|4]show quite a different 
picture. On the one hand, there is a strong correlation, indicated by distinct peaks, of 
angular position between hydrogen atoms in the same layer (Hi , H2 , and H3 positions 
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Figure 6. Absolute angular distribution of hydrogen atoms in the 2x2x2 
supercell. The result are obtained from a 39 ps MD run at 30 K. Positions of 
hydrogen atoms with labels H1,H2,- ■ -,1116 are shown in figure 1. 



are as in figure^. On the other hand, no visible correlations were found for the atoms 
in the next layer (H4, H5, and Hg). This means that the hydrogen atoms maintain 
relative orientations shown in figure [21 keeping their angular positions 120° apart — 
similar to the angular positions in the XGT/XLT structures mentioned in section [1] 
— but only if we consider atoms which belong to the same layer of the structure. 
On the contrary, the neighbouring layer positions appear only loosely coupled. These 
relations are maintained up to the highest temperature we have investigated (300 K), 
but broadening of the peaks observed at 244 K (upper panel of figure |4]) suggests that 
the peaks may disappear at higher temperatures, as one would expect. 

The above result supports, to some extent, the conclusions of the previous work on 
brucite crystal 13 [121 [HI HZl UHl US] j except that the hydrogen- hydrogen correlations 
do not appear for all atoms, but only for neighbouring atoms in the same layer. 
Furthermore, the lock-in of hydrogen positions to the lattice is not well confirmed by 
the above results and probably vanishes above relatively low temperature of T = 15 
K. These discrepancies motivated us to investigate the system further under less 
constraining conditions. 

The results of the calculation performed for the larger 2x2x2 unit cell at 
temperature T = 16 K are presented in figure [5j It is quite clear that the correlations 
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Figure 7. Relative angular distribution of hydrogen atoms in the 2x2x2 
supercell. The results are obtained from a 39 ps MD run at 30 K. Solid line 
shows the function Ta(}P) given by equation JSj fitted to the data. Positions of 
hydrogen atoms with labels H1,H2,- ■ -,1116 are shown in figure f . 



seen in the smaller unit cell are changed by the increase of the unit cell. We were unable 
to identify clearly any lock-in effect of the hydrogen positions with respect to the 
crystal lattice for the temperatures above 30 K (figure IH but below this temperature 
the signal arising from locked-in hydrogen positions is distinct and robust (figure [5]) . 
We have performed the calculations for several temperatures in the low temperature 
regime and for different starting conditions and obtained the same result in all cases. 
The relative angular correlations are still present when a larger unit cell is 
investigated, but they are also modified with respect to the ones found before with 
the smaller unit cell. The histograms seen in figure |6] show distributions of angular 
positions of hydrogen atoms relative to the crystal lattice and those shown in figure [7] 
— relative to an arbitrarily selected reference atom Hi, see figure [T] For the low 
temperature run (figure [5]) we can easily identify obvious non-uniformities in the 
distribution of absolute angular positions of most hydrogen atoms. The absence of 
distinct peaks for some atoms may be interpreted as a result of the relatively short 
simulation time (6 ps) and insufficient statistics due to the difficulty of performing 
adequate sampling of configuration space at low temperatures. We have verified that 
the overall picture remains in fact similar for different runs, but individual peaks 
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Figure 8. Relative angular distribution of hydrogen atoms in the 2x2x2 
supcrcell. The results are obtained from a 18 ps MD run at 150 K. Solid line 
shows the function Ta{}P) given by equation JSj fitted to the data. Positions of 
hydrogen atoms with labels H1,H2,- ■ -,1116 are shown in figure 1. 



disappear for different atoms. 

For higher temperatures, above 30-40 K, the absolute angular positions of 
hydrogen atoms are completely decoupled from the directions of crystal lattice 
(figureini) but we can still observe fairly strong correlations in relative angular positions 
of hydrogen atoms — these correlations are robust and persist even in the regime of 
rather high temperature. Figure [51 obtained for temperature T = 150 K, presents an 
essentially unchanged picture of the distribution of hydrogen atoms' relative positions. 
One can notice that weaker correlations with atoms from the first and from the last 
row seem to have disappeared. However, this may be due to the smaller statistics of 
the sample (18 ps vs. 39 ps) and higher noise in the data obtained for this case. 



J^.2. Angular distribution of hydrogen positions 

We have constructed a simple probabilistic model to describe measured angular 
distribution of hydrogen atoms. As a starting point it is reasonable to assume a 
Gaussian distribution around one location accompanied by a constant background 
filling for the rest of the histogram. In periodic angular coordinates the standard 
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Here 'd^ is an elliptic Jacobi function of the third kind. Thus, the angular distribution 
we introduce here reads: 



-F,(^)=AA 



(l_a)_L^3('^_Z,e-V)+« 



(3) 



where A/" is a normalizing factor. It consists of the correlated part given by equation 
^ with intensity (1 — a) and the background with intensity a. 

One can notice that in the case of 30 K (figure!?]) as well as for 150 K (figure[H]) the 
model for the angular distribution given by equation ([3]) fits the obtained numerical 
data remarkably well — it indicates the existence of a single potential minimum in 
relative motion of individual hydrogen atoms. This demonstrates that relative angular 
distribution refiects the correlated behavior of the hydrogen atoms, with a single 
optimal angular position for each hydrogen atom. This is obviously not the case 
in the absolute angular positions depicted in figure [SJ where a two-peaked distribution 
would be required to describe the data. 

We should also note that the atoms {Hg, Hio, Hii,Hi2} are characterized by 
strongly correlated movement with respect to the atom Hi. They are located in 
the same structure layer as the reference (Hi) atom. Also the position of a single 
out-of-plane atom H5, located directly above the reference atom Hi on the same 
three-fold axis, is strongly correlated with that of atom Hi. Moreover, it is clear that 
the directions of the peaks do not strongly follow 60° grid present in the smaller unit 
cell. The strongest signal indicates same/opposite direction arrangement of the direct 
neighbour (atom Hg) of the reference atom and its images in the neighbouring unit 
cells (H10-H12). However, if we look closer at the distributions presented in figure[71 we 
can notice small preference for image atoms (H2-H4) to cluster around 120° and 240° 
directions. This result does not offer strong support for the XGT/XLT arrangement 
hypothesis but the preferred directions are compatible with this type of arrangement. 

In light of these results we conclude that the angular correlations in the atomic 
distribution obtained for the small unit cell is likely just an artefact of the constraints 
imposed on the system by the -s/S x -\/3 x 1 unit cell. We also believe that it is still 
possible that the real physical system exhibits one of the symmetries proposed before 
O [7l [13l Uni [ITl HH] , but we emphasize that it does not show these symmetries for 
system size up to the 2x2x2 supercell. 

5. Summary and Conclusions 

We have performed extensive calculations to determine the structure and elastic 
properties of magnesium hydroxide. The static calculations confirm the main result 
of the previous research that static average structure with hydrogen atoms on the 
three-fold axis is absent in the brucite crystal. Further static calculations of the elastic 
constants of brucite also support conclusions from previous papers [HI [3 [131 HZl UHl HH] 
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suggesting existence of the superstructure, probably with the P3 symmetry, and XGT 
arrangement of the hydrogen atoms [13]. This, however, does not exclude all other 
possibilities which would be still consistent with the available experimental evidence. 

In the present paper we have investigated one of such possibilities — existence of 
hydrogen sublattice possibly almost decoupled from the rest of the crystal. The results 
from the MD calculations in the small supercell support this hypothesis. Furthermore, 
the calculations performed in the larger 2x2x2 cell revealed correlations between 
neighbouring hydrogen atoms at temperature between 30 K and 150 K. We have not 
found any lock-in behaviour of the hydrogen system above T = 30 K, while below 
this temperature we have identified strong correlations between angular positions of 
hydrogen atoms and orientation of the lattice vectors. This result indicates that the 
rotational (around the three-fold axis) degree of freedom of hydrogen system gets 
decoupled from the crystal lattice around temperature T = 30 K which may be 
considered as a characteristic energy scale for the onset of hydrogen-lattice coupling. 

The 2x2x2 supercell (shown in figure [T]) is the largest one available for practical 
use at the moment, and additional calculations of lattice dynamics performed for 
the 3x3x3 supercell indicate that the long-range interactions do not modify lattice 
dynamics in any significant way. Thus, we expect that the main qualitative conclusions 
on the crystal symmetry and hydrogen angular correlations either would not change, 
or would even strengthen in a larger system, when technology and algorithms will 
develop to the point that larger systems could be investigated. 

The results presented above demonstrate that the angular coordinates of the 
hydrogen atoms in the brucite crystal are likely decoupled from the crystal lattice 
above the estimated characteristic temperature T ~ 30 K. We have provided 
evidence that, in spite of this partial disorder in the hydrogen subsystem, angular 
correlations between positions of hydrogen atoms are robust and survive also in the 
high temperature regime T ~ 150 K. We believe that these results will stimulate 
further research and the issue of inter-hydrogen correlations in brucite and related 
crystals should be further investigated using theoretical as well as experimental 
methods. We suggest that nuclear magnetic resonance spectroscopy could be a useful 
experimental tool in this context. 
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